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We consider mesoscopic thermal transport between two bulk dielectrics joined by a narrow wire 
or weak mechanical link. In the "tunneling" regime where the phonon transmission probability 
through the link is small and the thermal conductance is much less than Tvk'^T/Gh, the thermal 
current is determined by a product of the local vibrational spectral densities of the two bodies. We 
derive an expression for the thermal current that is a thermal analog of the well-known formula for 
the electrical current through a tunneling barrier. 

PACS: 63.22.+m, 66.70.+f, 68.65.-k 
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I. INTRODUCTION 

Recently there has been considerable interest in the 
physics of phonons in mesoscopic and nanoscale systems. 
In a beautiful experiment, Schwab et al. W observed a 
quantization of the thermal conductance in freely sus- 
pended one-dimensional dielectric wires, corresponding 
to a conductance given by irk^T/GTi per transmitted vi- 
brational mode, analogous to the well-known electrical 
conductance quantization in units of e^ /2i:h per spin- 
resolved channel in onc-dimcnsional mesoscopic conduc- 
tors P]. The thermal conductance quantization can 
be understood f^-j^] by using a thermal analog of the 
Landauer-Biittiker formula |2| . 

In this paper we consider mesoscopic thermal transport 
through a wire or weak link in the limit where the phonon 
transmission probability through the wire is small (the 
so-called nonadiabatic regime) and the thermal conduc- 
tance is much less than nk^T/6h, the value correspond- 
ing to one fully propagating channel. The geometry we 
consider is shown schematically in Fig. nl 



The mechanical bridge may consist of one or more 
chemical bonds, or by a narrow wire of dielectric ma- 
terial, both of which we model by a harmonic spring of 
stiffness K [^. We obtain a general expression for the 
thermal current that can be regarded as a thermal ana- 
log of the well-known formula, derived by Schrieffer et al. 
[M,p|, for the electrical current through a tunneling bar- 
rier. Our result can also be interpreted as an application 
of the thermal Landauer formula in the weak-tunneling 
limit, with the energy-dependent phonon transmission 
probability calculated microscopically. 



II. THERMAL CURRENT IN THE 
WEAK-TUNNELING LIMIT 



We consider two macroscopic dielectrics, L and R, held 
at fixed temperatures Tl and Tr . The Hamiltonian of the 
isolated solids is Hq = Hi, + i^R, where 



Hr = 



7 ^^Jrt Qj„a/n, 



/ = L, R. 



(1) 
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FIG. 1. Mesoscopic wire connected nonadiabatically to 
bulk dielectrics. 



The al^j and a„/ are phonon creation and annihilation 
operators for the left and right sides. The vibrational 
modes of the isolated bodies are labeled by n and have 
energies ojin- 

The two bodies are connected by a harmonic spring 
with stiffness K, 



6H=\K{ul-ul 



(2) 



Here u| is the normal component of the displacement 
field u(r) at the surface of body / at the point of con- 
nection to the weak link, with the local surface normal 
taken to be in the z direction. Our mesoscopic weak-link 
model is illustrated in Fig. ||. 

The system is described by the Hamiltonian H = 
Hq + 5H . The thermal current is calculated by defin- 
ing a thermal current operator /th according to 



/th = 6)ti/R = i[H,HR]. 



(3) 



The expectation value of Ith is the energy per unit time 
flowing from the left to the right body. To leading order 
in SH (the weak-coupling limit) we obtain M 



2TrK^ 



Hh 



deeN^{e)Nn{e)[n^{e)-nn{e)], (4) 



an expression analogous to the formula derived by Schri- 
effer et al. MM for the electrical current through a tunnel- 
ing barrier. Here ni^{e) and rtR(e) are Bose distribution 
functions with temperatures Tl and Tr, and 



Ni{Lo) = -^lmDiiLo) 



(5) 



is the local spectral density, defined in terms of the 
Fourier transform of the retarded correlation functions 



Djit)^-t0{t){[u}{t),um]), 



(6) 



for the isolated macroscopic bodies L and R. The spectral 
density defined in dq) is different than the ordinary ther- 
modynamic DOS, even in a homogeneous system. The 
quantity (H), however, is the one relevant here. 



For example, a = 1 at the planar surface of a semi-infinite 
isotropic elastic continuum. 

A simple application of our theory to a nanostructure 
consisting of a cylindrical neck of Si material connecting 
two semi-infinite Si crystals has been given in Ref. M. 



III. DISCUSSION 

The classical theory of thermal conduction, based on 
the heat equation and on the concept of a local ther- 
mal conductivity, is entirely inapplicable to mesoscopic 
dielectrics. In a mesoscopic system, thermal resistance 
is caused by elastic scattering of phonons, whereas in 
an infinite, disorder-free crystal it is caused by anhar- 
monicity. In the non-mesoscopic regime (for example, at 
higher temperatures), the thermal conductivity of a nar- 
row bridge is determined by the bridge material's bulk 
thermal conductivity k, but in the mesoscopic regime it 
is determined by the mechanical properties of the bridge 
material, through the stiffness K. 
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FIG. 2. Our weak link model. Two macroscopic di- 
electrics are joined by a harmonic spring. 

Our result (|j) shows that the (linear) thermal conduc- 
tance of a weak mesoscopic link, defined by 



G 



th 



iim . 

Tl^Tr Tl-Tr' 



(7) 



will vary at low temperatures as a power-law in T. If 
we let a be the exponent characterizing the power-law 
phonon spectral density of the reservoirs at low energies, 
then 



Gth oc t2"+i. 



(8) 
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